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Abstract:  In  this  paper  we  present  a  new  combinatorial  problem,  called  minmax  multidimensional  knapsack  problem  (MKP), 
motivated  by  a  military  logistics  problem.  The  logistics  problem  is  a  two-period,  two-level,  chance-constrained  problem  with 
recourse.  We  show  that  the  MKP  is  NP-hard  and  develop  a  practically  efficient  combinatorial  algorithm  for  solving  it.  We  also 
show  that  under  some  reasonable  assumptions  regarding  the  operational  setting  of  the  logistics  problem,  the  chance-constrained 
optimization  problem  is  decomposable  into  a  series  of  MKPs  that  are  solved  separately.  ©  2007  Wiley  Periodicals,  Inc.  Naval  Research 
Logistics  54:  656-666,  2007 
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1.  INTRODUCTION 

We  present  a  new  combinatorial  model  called  minmax 
multidimensional  knapsack  problem  (MKP).  The  MKP  is 
motivated  by  a  military  logistics  optimization  problem  typ¬ 
ical  to  tactical-level  ground  operations  such  as  supply  and 
resupply  of  ammunition  to  an  artillery  battalion  in  a  2-day 
operation.  The  problem  is  a  two-period  stochastic  program¬ 
ming  problem  with  recourse.  However,  under  our  reasonable 
operational  assumptions,  it  is  shown  to  have  unique  features 
that  reduce  it  into  two  separate  and  sequential  (albeit  not  inde¬ 
pendent)  sets  of  MKPs.  We  show  that  the  MKP  is  NP-hard  and 
develop  a  new  practically  efficient  algorithm  for  solving  it. 

Stated  simply,  the  military  logistics  problem  is  to  find  min¬ 
imum  cost  inventories  that  satisfy  minimum  responsiveness 
requirements.  The  requirements  are  expressed  in  terms  of 
probabilities  for  satisfying  demands  in  a  two-period  operation 
of  a  military  unit  that  comprises  several  weapons.  While  some 
inventories  must  be  determined  before  the  operation,  others 
can  be  set  after  the  demands  in  the  first  period  are  realized — a 
situation  that  lends  itself  to  a  two-period,  chance-constrained, 
stochastic  programming  model  with  recourse. 

Stochastic  programming  models  with  recourse  have  been 
applied  in  supply-chain  and  related  problems,  mostly  in  the 
context  of  optimizing  expected  values.  Escudero  et  al.  [9] 
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utilize  scenario  modeling  for  production  and  capacity 
planning.  Their  multiperiod  model  minimizes  expected  costs, 
while  considering  several  recourse  alternatives.  Dempster 
et  al.  [7]  consider  a  multiperiod  supply-chain  scheduling 
problem  where  demands  and  costs  are  uncertain.  In  a  recent 
paper,  Cattani  et  al.  [4]  analyze  the  simultaneous  produc¬ 
tion  of  two  products:  a  market-specific  product  tailored  to  the 
needs  of  individual  regions,  and  a  global  product  that  could  be 
sold  in  many  regions.  They  consider  a  two-stage  model  with 
additional,  post-recourse,  uncertainty  in  the  second  stage,  and 
seek  to  maximize  an  expected  profit  function.  Cheung  and 
Powell  [6]  consider  the  class  of  multistage  dynamic  networks 
with  random  arc  capacities,  and  propose  a  successive  convex 
approximation  approach  for  the  expected  recourse  function, 
which  captures  the  future  effects  of  current  decisions  under 
uncertainty.  This  method  decomposes  the  network  in  each 
stage  into  sub-problems  for  which  expected  recourse  func¬ 
tions  are  easy  to  obtain  [5].  Recently,  Gupta  et  al.  [11]  used 
chance  constraints  programming  approach  coupled  with  two- 
stage  stochastic  programming  with  recourse  methodology  to 
construct  a  two-stage  supply  chain  plan  under  demand  uncer¬ 
tainty  with  continuous  random  variables.  The  authors  utilized 
linear  programming  duality  to  obtain  the  expectation  of  the 
recourse  function  associated  with  the  second  supply  chain 
stage  in  terms  of  the  first  stage  production  decisions. 

One-period  chance-constrained  (OPCC)  problems  of  the 
form  min{cx  :  P[Tx  >  q]  >  p.  Ax  >  b,  x  >  0}  are 
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well  studied  for  continuous  random  variables  q.  However, 
there  are  only  a  few  papers  investigating  the  discrete  random 
variables  case,  e.g.  [3,8, 16],  In  [16]  the  author  considers  an 
integer  version  of  the  OPCC  problem  (IP-OPCC)  with  more 
general  objective  functions,  and  uses  methods  of  disjunctive 
programming  to  approximate  the  convex  hull  of  the  feasi¬ 
ble  region.  For  some  special  cases  a  full  description  of  the 
convex  hull  is  given.  In  [8]  the  concept  of  p-efficient  points 
(PEP)  is  introduced  and  used  for  deriving  lower  and  upper 
bounds  for  the  OPCC  problem  with  discrete  probability  dis¬ 
tributions.  Let  F  denote  the  joint  distribution  function  of  q, 
then  for  0  <  P  <  1,  the  PEP  are  defined  as  the  minimal 
realization  q  points  (in  the  sense  of  the  order  <)  for  which 
F(q)  >  p.  In  [2]  a  branch  and  bound  algorithm  that  utilizes 
PEP  is  presented.  In  [  15]  a  generalization  of  the  IP-OPCC  is 
analyzed  where  new  valid  inequalities  related  to  the  prece¬ 
dence  constrained  knapsack  polyhedra  are  developed  and  are 
used  by  a  general  iterative  algorithm  to  solve  the  general¬ 
ized  problem.  In  these  methods  an  optimal  solution  can  be 
obtained  by  exploring  the  set  of  PEP.  However,  since  the  set  of 
PEP  can  be  extremely  large,  the  optimization  procedure  may 
encounter  computational  difficulties.  Nevertheless,  “good” 
solutions  can  be  derived  by  looking  at  a  subset  of  PEP. 

In  this  paper  we  use  a  different  IP  formulation  to  model 
the  OPCC  problem  with  discrete  probability  distributions, 
and  observe  that  the  resulting  IP  problem  is  a  special  case  of 
a  combinatorial  problem  termed  minmax  multidimensional 
knapsack  problem  (MKP).  Given  a  set  A  —  {as  :  s  e  S’] 
of  n -dimensional  vectors  a'  of  real  numbers,  each  associ¬ 
ated  with  a  weight  ps ,  the  MKP  consists  of  selecting  a  subset 
A'  C  A  such  that  ^aSgA,  ps,  the  weight  of  A',  exceeds  a 
certain  threshold,  and  the  objective  is  to  minimize  the  sum  of 
the  componentwise  maxima  of  the  vectors  in  A'. 

The  one-dimensional  MKP  is  related  to  the  min  knapsack 
problem  [1],  where  one  looks  for  a  subset  of  items  such  that 
the  sum  of  their  weights  exceeds  a  given  constant,  and  the  sum 
of  their  values  is  minimized.  While  both  problems  have  the 
same  cover  type  constraint,  they  differ  in  the  objective  func¬ 
tion;  in  the  MKP  we  look  for  a  minmax  solution  while  in 
the  min  knapsack  problem  we  seek  a  minimum  sum  solu¬ 
tion.  The  different  objective  functions  affect  the  complexity 
of  the  problems;  while  the  min  knapsack  problem  is  known  to 
be  NP-hard  (see  for  example  [1])  the  one-dimensional  MKP 
can  be  solved  by  a  simple  greedy  polynomial  time  algorithm. 
Another  remotely  related  problem  is  the  NP-hard  multidi¬ 
mensional  knapsack  problem  which  also  consists  of  selecting 
A',  a  subset  of  a  given  set  of  items,  such  that  the  total  value 
of  A'  is  maximized  while  a  set  of  knapsack  constraints  are 
to  be  satisfied  (e.g.  [14]  and  [3]).  We  note  that  the  multidi¬ 
mensional  knapsack  problem  differs  from  the  MKP  in  both 
the  constraints  and  the  objective  function.  The  constraint  in 
the  MKP  is  of  a  covering  type  while  the  constraints  of  the 
multidimensional  knapsack  problem  are  of  a  packing  type. 


Also  the  objective  function  of  the  multidimensional  knap¬ 
sack  problem  is  to  maximize  a  total  value  while  in  the  MKP 
it  is  a  minmax  objective  function.  To  the  best  of  our  knowl¬ 
edge  this  is  the  first  time  the  MKP  is  defined  and  solved. 
We  believe  that  the  MKP  and  its  algorithm  are  interesting  by 
themselves  and  moreover,  are  useful  for  solving  a  class  of 
chance-constrained  problems. 

The  rest  of  the  paper  is  organized  as  follows.  In  Section  2 
we  motivate  the  modeling  effort  and  introduce  some  nota¬ 
tion.  In  Section  3  we  introduce  the  single-period  chance- 
constrained  optimization  problem,  and  in  Section  4  we 
describe  the  MKP  and  develop  an  algorithm  for  solving  it. 
We  show  that  the  single-period  chance-constrained  prob¬ 
lem  is  an  MKP.  In  Section  5  we  formulate  the  two-period 
chance-constrained  optimization  problem  and  show  that  this 
formulation  is  equivalent  to  an  IP  model,  which  is  decom¬ 
posed  into  a  series  of  MKPs.  Summary  and  conclusions  are 
given  in  Section  6.  The  Appendix  contains  proofs  for  the 
decomposition  property. 


2.  MOTIVATION  AND  NOTATION 

Our  model  is  motivated  by  a  military  logistics  problem 
typical  to  ground  operations.  Consider  a  battalion  that  com¬ 
prises  several  weapons  (e.g.,  artillery  pieces).  Each  weapon 
consumes  ammunition  from  a  designated  attached  stockpile. 
There  is  also  a  considerable  amount  of  ammunition  on-board 
the  weapon  itself  [13];  however  this  inventory  is  an  emer¬ 
gency  safety  stock,  to  be  used  only  if  the  stockpile  is  empty. 
The  total  amount  of  ammunition  in  the  stockpile  and  on  board 
the  weapon  is  assumed  to  be  enough  to  satisfy  any  foresee¬ 
able  demand  but  the  objective  is  to  tap  the  on-board  inventory 
as  little  as  possible  and  to  rely  only  on  the  stockpile  [10]. 
An  additional  ammunition  depot  is  attached  to  the  battalion 
headquarters  at  the  rear  of  the  combat  zone. 

We  consider  a  2-day  operation.  During  the  first  day  a 
weapon  can  use  ammunition  only  from  its  attached  stock¬ 
pile  or,  if  needed,  from  its  own  safety  stock.  It  cannot  rely 
on  logistical  support  from  other  sources  because  movement 
in  the  combat  zone  during  the  combat  operation  is  risky.  At 
the  end  of  the  first  day,  after  demands  have  been  observed, 
depleted  stockpiles  and  possibly  reduced  safety  stocks  are 
replenished  from  the  depot  or  from  other  weapons’  stock¬ 
piles  that  transship  surplus  ammunition.  The  replenishment 
process  is  completed  before  the  second  day  of  operation. 
The  replenished  amount  of  each  weapon  must  cover  pos¬ 
sible  expenditure  from  the  weapon’s  safety  stock  during  the 
first  day  and  the  (yet  unknown)  demand  in  the  second  day. 
Similar  inventory  control  situations  may  occur  in  the  retail 
industry  too. 

The  weapons  are  called  henceforth  demand  points  (DP). 
We  require  that  the  supplies  allocated  to  the  DPs  on  each  one 
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of  the  two  days  satisfy  the  demand  during  that  day  with  a 
given  minimum  probability.  Because  of  the  inherent  scarcity 
of  relevant  demand  data  regarding  military  operations  [12], 
demand  distributions  are  generated  based  on  expert  inputs 
in  the  form  of  combat  scenarios.  In  the  one-period  prob¬ 
lem,  let  ZD  =  (D\,  £>2, . . . ,  Dn)  denote  the  vector  of  random 
variables  that  represents  the  demands  at  DP,,  i  —  1 ,...,«, 
and  let  d  —  (d\ , . . . ,  dn)  denote  a  realization  of  this  vector, 
which  is  called  a  demand  scenario.  We  assume  a  finite  num¬ 
ber  of  demand  scenarios  where  S  denotes  their  index  set.  Let 
Pr[D  —  ds]  =  ps  >  0,j  e  S,  where  P*  ~  1-  In 
the  two-period  problem  we  have  two  vectors  of  random  vari¬ 
ables,  and  accordingly,  two  index  sets  of  demand  scenarios, 
Sl  and  S2. 

Let  Xj  denote  the  amount  of  supply  initially  allocated  to 
D  Pj ,  and  let  Y  denote  the  amount  of  supply  initially  deployed 
in  the  depot.  The  parameters  CXj  and  Cy  are  the  cost  of  a  unit 
of  supply  at  DP,  and  the  depot,  respectively.  These  costs  are 
incurred  at  the  first  period.  There  are  no  additional  costs  in  the 
problem.  Military  DPs  are  similar  units  (e.g.,  artillery  guns) 
and  therefore  it  is  reasonable  to  assume  that  CXi  =  C x  for  all 

1  =  1 ,...,«. 

The  decisions  made  at  the  beginning  of  day  1  are  with 
respect  to  (a)  the  amount  of  supply  x;  (stockpile)  to  be  allo¬ 
cated  to  each  DP,  and  (b)  the  total  amount  of  supply  Y  to 
be  kept  in  the  depot  at  the  beginning  of  period  1,  to  be  used, 
if  necessary,  at  the  second  period.  The  recourse  variables 
are  the  shipments  that  take  place  at  the  beginning  of  day 

2  from  and  to  the  depot.  For  a  given  pair  of  distributions  of 
interrelated  discrete  demand  random  variables — one  for  each 
day — the  objective  is  to  find  a  minimum  cost  inventory  policy 
such  that  demands  are  satisfied  in  both  days  with  probabili¬ 
ties  that  exceed  certain  thresholds.  From  now  on  we  use  the 
terms  “day”  and  “period”  interchangeably.  We  begin  with  the 
single-period  problem. 


3.  THE  SINGLE  PERIOD  PROBLEM 

The  problem  is  to  minimize  the  total  amount  of  supply  such 
that  a  certain  level  of  logistics  responsiveness  is  attained.  The 
logistic  responsiveness  is  measured  by  the  probability  that  all 
demands  are  satisfied,  that  is,  none  of  the  weapons  has  to  use 
ammunition  from  its  safety  stock.  Formally, 


n 


min  xi 

7  =  1 

(1) 

s.t. 

Pr[D  <  x]  >  Q, 

(2) 

x  >  0,  x  eZn. 

(3) 

where  x  =  (x\,.. .  ,xn). 
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Constraint  (2)  is  the  logistics  responsiveness  chance 
constraint  and  Q  is  the  probability  threshold  set  by  the  com¬ 
mander.  Since  we  assume  that  the  demand  scenarios  are  dis¬ 
crete  random  variables,  problems  (l)-(3)  may  be  formulated 
as  the  following  IP  problem: 

n 

min  ^2  Xi  (4) 

1  =  1 
s.t. 

xi  —  d-Ss  >  0,  i  =  n,  s  e  S,  (5) 

^2  P'Ss  -  Q’  (6) 

seS 

Xi  >  0,  X{  e  Z,  Ss  e  [0, 1},  i  =  1, . . .  ,n,  s  e  S.  (7) 

The  binary  variables  Ss,  s  e  S,  indicate  whether  a  cer¬ 
tain  scenario  has  been  selected  to  be  satisfied.  Constraint  (6) 
guarantees  that  the  probability  of  unsatisfied  demand  does 
not  exceed  the  operationally  set  threshold  1  —  Q. 

Problems  (4)— (7)  is  an  MKP,  a  newly  defined  combinato¬ 
rial  problem  which  is  formulated  next. 


4.  THE  MKP 

Let  as  —  (t«j, . . . ,  asn)  be  an  n -dimensional  real  vector,  and 
let  A  =  {as  :  s  e  S}  be  a  set  of  |S|  such  vectors.  Each  vector 
as  in  A  has  weight  (value)  ps  >  0.  Let  A!  C  A  be  a  subset 
of  A.  Define 

n 

G(A')  —  max  a?, 

z '  aseA' 

1  =  1 

where  maxa>e/v  a]  is  the  largest  ith  component  among  the 
vectors  in  A! . 

A  subset  A!  C  A  is  said  to  be  Q-feasible  for  a  given 
parameter  Q  if  ps  >  Q.  The  minmax  multidimen¬ 

sional  knapsack  problem  (MKP)  is  to  find  a  <2-feasible  subset 
A*  C  A  such  that  G(A*)  is  minimal.  The  subset  A*  is  called 
an  optimal  minmax  subset. 

EXAMPLE!:  Let  A  =  {(1, 5),  (5, 3),  (3, 7),  (10, 1),  (7, 7)}, 
Q  —  0.8  and  ps  —  0.2  for  all  s.  Then  any  subset  of  A 
containing  at  least  4  elements  is  feasible,  and  an  optimal  min¬ 
max  subset  of  A  is  A*  =  {(1, 5),  (5, 3),  (3, 7),  (7, 7)},  with 
G(A*)  =  7  +  7=  14. 

It  is  easily  seen  that  problem  (4)— (7)  is  actually  an  MKP, 
where  as  is  a  demand  scenario  ds,  s  e  S,  and  the  objective 
is  to  find  a  0-feasible  subset  of  these  scenarios  such  that 
the  sum  of  the  componentwise  maxima  of  the  corresponding 
demand  vectors  is  minimized. 
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Next  we  develop  a  simple  combinatorial  algorithm, 
termed  the  MKP-Algorithm,  for  solving  the  MKP.  The  case 
n  =  1  is  trivial.  A  greedy  algorithm  solves  the  problem 
by  simply  arranging  the  elements  in  A  in  a  nondecreas¬ 
ing  order  .y(l), . . . ,  s(|S|).  The  optimal  subset  is  A*  = 
as(2\  . . .  ,as V  '},  where  j*  is  the  smallest  index  such 
that  ps(j)  >  Q. 

The  case  n  >  1  is  more  complex.  The  algorithm  com¬ 
prises  two  stages:  a  preliminary  stage  and  a  main  stage.  In 
the  preliminary  stage  we  discard  from  A  all  vectors  that  are 
evidently  nonoptimal.  In  the  main  stage  the  remaining  sub¬ 
set  is  "pruned”  and  an  optimal  minmax  subset  is  obtained. 
That  is,  the  algorithm  finds  and  then  deletes,  a  complement 
of  a  (^-feasible  subset,  to  obtain  a  required  optimal  minmax 
subset. 


4.1.  Preliminary  Stage 

Let  «s'um  =  Yl'i=iai’s  e  5,  and,  without  loss  of  generality, 
assume  that  a*  <  a2  <  ...  <  aJf'.Let  j*  be  the  smallest 
index  such  that  ^=1  ps  >  Q-  Define  the  set  B  as  follows: 

B  =  J«s  e  A  :  assum<J2  max  W>  •  •■>«/} 


Note  that  {a 1 , . . . ,  aJ‘  }  is  a  subset  of  B.  For  an  empty  set  we 
let  j*  =  0. 


PROPOSITION  4.1:  If  as  e  A  \  B  for  some  ,v  e  S  then 
as  A*,  where  A*  is  any  optimal  minmax  subset. 

PROOF:  Suppose  as'  e  A  \  B  and  as'  e  A*  for  some 
optimal  minmax  subset  A*.  Then, 


G(A*)  =  maxaf  >  «sum- 
z — '  aseA* 
i= 1 

Since  as'  e  A  \  B,  it  follows  that 

n 

flsum  >  max  {ai  >■■■■>  ai  }  =  G({a\...,aJ‘}), 

1=1 

where  {a1, . . . ,  a7*}  is  a  feasible  subset.  Thus,  A*  cannot  be 
an  optimal  minimax  subset,  in  contradiction.  □ 

We  conclude  that  an  optimal  minmax  subset  of  A  must  be 
a  subset  of  B. 


4.2.  Main  Stage 

We  search  for  an  optimal  minmax  subset  by  considering 
its  possible  complement  subsets  in  B.  In  order  to  reduce  the 
number  of  subsets  examined,  we  generate  T,  a  family  of 
subsets,  in  which  a  complement  of  an  optimal  minmax  sub¬ 
set  is  guaranteed  to  be  included.  We  describe  the  construction 
of  'P  more  precisely  below. 

For  each  i,i  =  1, ...» n,  we  arrange  the  vectors  of  B  in  a 
nonincreasing  sequence  B,  according  to  their  ith  component. 
That  is, 

as‘(r>  >  cF,(2)  >-•••>-  aSi(IN)  if  anc|  oniy  jf  as.^1) 


and  B,  is  the  sequence  as,(1\  . . . 

Then,  we  define  T  to  be  the  family  of  all  possible  sub¬ 
sets  ijfk,k  —  1, . . . ,  I'Ll,  of  B  that  satisfy  the  following  two 
conditions: 

(1)  t/r*  is  of  the  form 

fk=  U 

I€{1, 

where  0  <  m,  <  | B\.  Thus,  each  ^  is  a  pre¬ 
fix,  that  is,  a  union  of  truncated  sequences  from  /!,, 
i  =  1, . . . ,  n.  Note  that =  0  means  that  no  prefix 
is  taken  from  the  i  th  order. 

(2)  xjfk  is  a  maximal  cardinality  set  that  satisfies 

p(A  \B)+  Y2  P>  -  1  ~  G’ 

a‘£\lsk 


where  p(A\B)  =  X!ase(A\B)  Ps  •  That  is,  the  comple¬ 
ment  in  A  of  the  union  of  the  set  of  deleted  vectors 
(which  are  evidently  nonoptimal)  and  x[rk,  is  a  Q- 
feasible  subset  of  A.  Adding  one  more  vector  from 
B  to  x/fk  will  violate  the  inequality  and  therefore  will 
render  the  complement  infeasible. 

Now,  for  each  x(rk  e  'P  define 

n 

M{fk)  =  V  max  a-  —  G(B  \  fk) 


and  let 


k  e  argmin{M(i/fi)  :  fa  e  'P}. 

The  set  B  \  \jj-k  is  an  optimal  minmax  subset  of  A.  Before 
we  prove  the  validity  of  the  MKP-Algorithm,  consider 
Example  2. 
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Table  1.  Scenarios. 


5 

a\ 

a2 

1 

3 

5 

2 

6 

1 

3 

8 

9 

4 

6 

5 

5 

3 

5 

6 

4 

9 

7 

2 

7 

8 

7 

2 

9 

10 

10 

10 

4 

3 

EXAMPLE  2:  Let  S  =  {1, ... ,  10},  n  =  2,  Q  =  0.6,  and 
ps  =  0.1  for  all  s  e  S  (see  Table  1). 


First,  we  rank  the  vectors  in  a  nondecreasing  order  of  the 
sum  of  their  components  (see  Table  2). 

j*  =  6;  Ml)  =  2,s  (2)  =  10,5(3)  =  1, 

5(4)  =  5,5(5)  =  7,5(6)  =  8; 

{7  10  1  5  7  R 1 

a\,ax  9a[, 

+  max  {ala^,a\,al,al,al}  =  7  +  7  =  14. 

Proposition  4.1  shows  that  a3  and  a 9  can  be  deleted  from 
further  consideration  since  the  sums  of  their  components  are 
17  and  20,  respectively,  which  are  greater  than  14.  Thus, 
B  —  {a1, a2, a4, a5, a6, a7, a8, a10}. 

Next  we  generate  two  rankings  of  B,  according  to  each 
component  of  as .  Each  ranking  is  in  a  nonincreasing  order  of 
the  corresponding  component  (see  Table  3). 

Now,  *  =  {{(7, 2),  (6, 1)},  {(7, 2),  (4, 9)},  {(4, 9),  (2, 7)}}, 
and  Af({(7,2),(6, 1)})  =  6  +  9  =  15,  M({(7,2),  (4,9)})  = 
6  +  7  =  13  and  M({( 4, 9),  (2, 7)})  =  7  +  5=12. 

We  conclude  that  the  optimal  minmax  subset  A*  comprises 
the  vectors  with  the  indices  1, 2, 4, 5,  8, 10,  and  G(A*)  =  12. 
Note  that  the  optimal  minmax  subset  in  this  example  is 
unique.  This  is  not  necessarily  true  in  general. 

Recall  from  Proposition  4. 1  that  any  optimal  minmax  sub¬ 
set  of  A  is  contained  in  B.  The  next  proposition  shows  the 
existence  of  an  optimal  subset  such  that  its  complement  in  B 
lies  in  + . 

PROPOSITION  4.2:  There  exists  an  optimal  minmax 
subset  A*  c  B  such  that  B  \  A*  e  'P. 

PROOF:  Let  W  C  B  be  the  complement  of  an  optimal 
minmax  subset  in  B,  and  suppose  that  W  ^  +.  Define, 

E(W)  =  {a1  e  W  :  Vi  there  exists 

as‘  e  (B  \W)  and  af  >  a- }. 


That  is,  if  as  e  E(W),  then  for  each  B, ,  at  least  one  of  the 
components  of  the  truncated  sequence  (as,31),as'(2\  . . . 
is  not  contained  in  W.  Let  VV  be  the  set  of  all  such  W  and  let 
W  be  such  that 

\E(W)\  =  min{|£(W)|  :  W  e  WJ. 

If  |£(W)|  =  0  we  are  done.  Otherwise,  there  is  a  vector 
as  e  E(W).  According  to  the  definition  of  E(W),  for  each 
i  —  1, ....  n  there  exists  aSi  e  B  \  W  such  that  af  >  But 
then, 

n 

Y'  max  {«;  : 

1  =  1 

n 

=  Y. max  ( { aj  '■  e  B  \  W}  U  {«/}), 

i= 1 

which  implies  that  E(W)  is  not  minimal,  in  contradiction. 
Therefore,  there  exists  an  optimal  minmax  subset  whose 
complement  belongs  to  'P .  □ 

Thus,  the  algorithm  shown  above  solves  the  MKP.  Recall 
that  as  defined  in  [8],  for  any  distribution  function  E  otV/,  and 
any  0  <  p  <  1,  the  /^-efficient  points  (PEP),  are  defined  as 
the  minimal  realization  q  points  (in  the  sense  of  the  order  <) 
for  which  F(q)  >  p.  Notice  that  any  PEP  is  p-feasible  but  not 
vice  versa.  The  MKP  can  be  solved  by  the  more  general  and 
powerful  methods  which  are  based  on  polyhedral  techniques 
as  described  previously  in  the  Introduction  and  presented  in 
[2,  8,  15, 16].  However,  our  MKP-Algorithm,  which  exam¬ 
ines  points  in  B — possibly  some  non-PEP — is  conceptually 
simpler  and  easier  to  use  than  those  methods. 


4.3.  Complexity  and  Computational  Analysis 

THEOREM  4.3:  MKP  is  NP-complete  and  the  complexity 
of  the  MKP-Algorithm  is  0(tik),  where 


k  =  min 


~  1-6" 
min,  ps 


Table  2.  Rank  order  of  the  sums. 


s 

a2 

<um 

2 

6 

1 

7 

10 

4 

3 

7 

1 

3 

5 

8 

5 

3 

5 

8 

7 

2 

7 

9 

8 

7 

2 

9 

4 

6 

5 

11 

6 

4 

9 

13 

3 

8 

9 

17 

9 

10 

10 

20 
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PROOF:  We  prove  the  NP-completeness  by  showing  a 
simple  reduction  to  the  knapsack  problem.  Recall  that  the 
knapsack  problem  is  defined  as  follows.  Given  non-negative 
numbers  c,-,u>,-,  i  =  1 ,...,«,  and  b,  find  a  subset  T  C 
{1 such  that  r  u>j  <  b  and  r  c is  maximized. 
We  show  that  the  knapsack  problem  polynomially  transforms 
to  MKP  using  the  following  reduction.  Let  |S|  =  n  and  let  A 
be  the  set  of  n  vectors  of  dimension  n,  where  the  /  th  vector  is 
set  to  be  c,-  x  e,  with  e,  the  i  th  unit  vector.  Let  p,  =  w,  and 
Q  —  Y^'i=  i  wi  ~  b.  Then,  if  A*  c  A  is  an  optimal  ^-feasible 
subset  then  A  \  A*  is  an  optimal  solution  for  the  knapsack 
problem,  and  vice  versa.  We  turn  now  to  show  the  complexity 
of  the  MKP- Algorithm. 

Summing  up  the  n  components  for  each  one  of  the  vec¬ 
tors  in  A  takes  0(n\S\)  computations,  arranging  the  sums 
in  a  nonincreasing  sequences  takes  0(|S'|  log  |S|)  steps  and 
constructing  the  set  B  takes  (9  ( |  .S’  | )  steps.  Thus,  the  prelim¬ 
inary  stage  of  the  algorithm  takes  0(n\S\)  +  <9(|S|  log  |5|) 
steps.  The  main  stage  of  the  algorithm  requires  arranging  n 
sequences,  an  operation  that  requires  0(n\S\  log  |5’|)  com¬ 
putations.  Let  k  be  an  upper  bound  on  the  size  of  the  largest 
set  in  +  .  Observe  that  k  is  the  minimum  between  the  num¬ 
ber  IS |  of  possible  vectors  in  A,  and  the  maximal  number  of 
vectors  that  can  “squeeze”  in  any  i/q,.  This  maximal  number 
is  determined  by  the  smallest  value  among  { ps  :  s  e  S|  and 
the  threshold  Q,  and  is  given  by  f  ] .  Finally,  examining 
all  possible  subsets  in  +  takes  0(nk).  This  can  be  done  by 
considering  vectors  of  length  of  at  most  k ,  in  which  the  j  th 
component  is  a  number  1  <  i  <  n  indicating  according  to 
which  prefix  the  vector  (scenario)  was  chosen  to  be  included 
in  4^ .  Also,  after  each  step,  the  set  of  possible  vectors  (scenar¬ 
ios)  in  i/q,  is  updated  by  removing  the  already  chosen  vectors 
(scenarios).  Thus,  since  0(nk)  >  0(n\S\  log  |S|),  the  com¬ 
plexity  of  the  MKP- Algorithm  is  (9(«|S|)  +  0(|S'|  log  |  S  | )  + 
0(n\S\  log  |5|)  +  0(nk)  =  0(iik).  □ 

Note  that  examining  all  possible  subsets  of  +  can  also  be 
done  by  checking  all  prefixes  of  length  0  —  k  for  each  i .  This 
takes  a  time  of  at  most  0(kn)-  However,  since  in  practical 
military  (and  many  commercial)  logistics  problems  A:  is  a  rela¬ 
tively  small  integer,  we  have  chosen  to  present  the  complexity 


Table  3.  Rank  order  according  to  the  first  component. 
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a2 
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3 
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Table  4.  Average  running 
Algorithm  and  the  IP  code. 

times 

(in  seconds)  of  the 

MKP- 

MKP 

IP 

\s\ 

Average 

S.D. 

Average 

S.D. 

40 

0.07 

0.06 

0.62 

0.12 

60 

1.44 

2.66 

3.88 

2.09 

100 

71.41 

76.79 

375.67 

392.10 

in  terms  of  0(nk)  rather  than  0(k")-  The  probability  thresh¬ 
old  Q  is  usually  larger  than  0.9,  and  relevant  scenarios  have 
typically  probability  that  is  at  least  0.01.  Therefore,  in  such 
situations  k  <  10,  which  implies  the  high  efficiency  of  our 
algorithm  for  any  practical  purposes.  Note,  as  well,  that  for  a 
fixed  A:  or  a  fixed  n,  our  MKP- Algorithm  runs  in  polynomial 
time. 

Next,  the  MKP- Algorithm  is  compared  to  a  general  IP 
code  (CPLEX  8.0)  with  respect  to  running  time.  The  com¬ 
parison  has  been  executed  on  Intel  Pentium  4,  2Ghz  CPU, 
512Mb  RAM,  run  under  Windows  XP.  The  MKP- Algorithm 
has  been  implemented  in  Microsoft  Visual  C++  6.0. 

We  fix  the  number  of  DPs  at  n  =  15,  which  is  approxi¬ 
mately  the  number  of  battalions  in  a  division,  the  number  of 
batteries  in  an  artillery  regiment,  and  the  number  of  artillery 
pieces  in  a  battalion.  The  number  of  scenarios  ranges  between 
|S|  =  40  and  |S|  =  100.  The  scenarios  are  uniformly 
distributed  (all  scenarios  are  equally  probable)  and  the  prob¬ 
ability  threshold  is  0.9.  For  each  value  of  |5|,  100  randomly 
generated  problems  were  solved  by  the  MKP-Algorithm 
and  the  general  IP  code.  Table  4  presents  the  average  run¬ 
ning  times  (in  seconds)  and  the  corresponding  standard 
deviations. 

For  this  range  of  data  the  MKP-Algorithm  clearly  outper¬ 
forms  the  general  IP  code.  This  result  is  reversed  if  n  gets 
larger.  For  |S|  —  n  =  50,  the  general  IP  code  is  10  times 
faster  than  the  MKP-Algorithm.  However,  in  practical  mili¬ 
tary  logistics  problems,  in  which  n  is  relatively  small  (e.g.  an 
artillery  battalion  comprises  16  pieces),  the  MKP-Algorithm 
is  more  efficient. 


5.  THE  TWO-PERIOD  MODEL 

Recall  that  x,-  is  the  amount  of  supply  initially  allocated  to 
I)  P, ,  i  =  1 and  Y  is  the  amount  of  supply  deployed 
initially  in  the  depot.  These  are  the  first-stage  decision  vari¬ 
ables.  The  second  stage  (recourse)  variables  are  shipments 
between  the  depot  and  the  DPs  that  take  place  after  the  first 
day  demands  are  realized.  Suppose  that  scenario  s  e  S]  has 
been  realized  in  the  first  day.  For  that  particular  scenario  .v ,  and 
the  set  of  demand  scenarios  for  day  2,  with  their  correspond¬ 
ing  probabilities,  let  yf,  i  —  denote  a  feasible  flow 
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of  supply  between  the  depot  and  D  P;  such  that  the  demand 
in  day  2  is  satisfied  with  probability  not  smaller  than  a  given 
threshold.  If  y?  >  0,  then  the  depot  sends  out  supply  to  DPj, 
and  if  y.  <  0,  the  depot  receives  surplus  supply  back  from 
DP,.  The  sum  y5  =  ^>~"=|  y?  denotes  the  total  net  flow  of 
supply  between  the  depot  and  the  DPs  at  the  end  of  the  first 
day,  given  scenario  s  has  been  realized.  This  amount  is  to 
be  shipped  out  before  the  beginning  of  day  2,  to  satisfy  the 
demand  in  the  second  day.  If  ys  <  0,  then  there  is  a  back 
flow  of  supply  from  the  DPs  to  the  depot  at  the  end  of  the 
first  day,  which  means  that  the  supplies  left  in  the  DPs  after 
the  first  day  are  more  than  enough  for  satisfying  the  demand 
in  day  2.  Let  y  —  max^i})^}.  If  y  is  positive,  then  it  is  the 
amount  of  supply  that  is  stored  in  the  depot  at  the  beginning 
of  day  1,  that  is,L  =  max{0,  y}. 

Consider  the  following  example.  Suppose  n  —  2  and 
there  are  three  equally  likely  demand  scenarios  for  day  1 : 
(100,0),  (0, 150),  and  (30,40)  for  s  —  1,2,3,  respectively. 
These  three  demand  scenarios  must  be  satisfied  with  certainty 
(<2i  =  1).  Suppose  that  for  day  2  there  are  two  equally  likely 
scenarios:  (80,0)  and  (0,70)  that  must  be  satisfied  with  cer¬ 
tainty.  Clearly,  xi  =  100  and  xy  =  150  is  a  feasible  (and 
minimum)  solution  for  the  first  day.  Suppose  scenario  1  has 
occurred,  then  yj  =  80  and  y\  =  —80  is  a  feasible  flow 
of  supply  between  the  depot  and  the  two  DPs  at  the  end  of 
day  1 .  In  fact,  these  flows  are  also  minimal  in  the  sense  that  the 
depot  must  ship  at  least  80  units  to  DP\ ,  and  it  cannot  receive 
more  than  80  units  from  DPs.  Here  y1  =0.  If  scenario  2  has 
occurred  then  y\  —  —20  and  y\  —  70  with  y2  =  50  are  min¬ 
imum  feasible  flows,  which  imply  that  the  depot  must  have  at 
least  50  units  in  order  to  satisfy  the  required  demand  in  day  2. 
If  scenario  3  has  occurred,  then  y3  =  10  and  y|  =  —40 
with  y3  =  —30.  In  this  case  ( s  —  3)  the  net  flow  is  neg¬ 
ative  and  therefore  there  is  no  need  for  supplies  stored  in 
the  depot.  In  this  example,  y  =  max{0, 50,  —30}  =  50  and 
Y  —  max{0, 50}  =  50. 

For  simplicity  of  exposition,  as  long  as  it  does  not  cause 
any  confusion,  we  have  chosen  to  indicate  by  the  index  s 
the  variables  and  the  data  that  correspond  to  the  first  day 
and  by  the  index  t  those  that  correspond  to  the  second  day. 
Accordingly,  pst  is  the  conditional  probability  of  scenario 
t  e  S2  given  that  scenario  .y  e  S 1  has  been  realized  in 
day  1.  The  two-period  problem  is  described  as  follows:  In 
period  1,  before  observing  the  demands  in  day  1,  we  deter¬ 
mine  the  inventories  x,,  i  =  1, in  the  DPs  and  the 
inventory  Y  in  the  depot.  At  the  end  of  day  1,  after  observ¬ 
ing  the  demand  scenario  s  that  has  been  realized  in  day 
1,  we  determine  the  shipments  y?  to  and  from  the  depot 
to  be  carried  out  before  day  2.  Excess  supply  at  a  cer¬ 
tain  DP  is  shipped  back  to  the  depot,  and  shortage  in  sup¬ 
ply  at  a  ceratin  DP  is  replenished  from  the  depot.  The  net 
flow  of  these  shipments  determine  the  inventory  Y  in  the 
depot,  which  operates  also  as  a  clearing  house  for  the  DPs. 
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Recall  that  there  are  no  shipment  costs.  The  initial  alloca¬ 
tion  to  the  DPs  Xi,  i  —  1 must  satisfy  the  demands 
in  day  1  with  probability  not  smaller  that  Qi  -  see  con¬ 
straints  (9)-(ll)  below.  For  each  realized  scenario  in  day 
1,  including  a  scenario  s'  for  which  df  >  Xj  for  some  i  (in 
which  case  the  safety  stock  in  DP,  is  used),  if  such  exists, 
we  require  that  (a)  all  on-board  safety  stocks  are  replen¬ 
ished  to  their  full  initial  capacity,  and  (b)  the  demands  in 
day  2  are  satisfied  with  probability  not  smaller  than  Qy. 
These  requirements  are  manifested  in  constraints  (12)— (14) 
below.  The  two-stage  optimization  model  is  formulated  as 
follows: 


min  CxX  +  CyY 

(8) 

s.t. 

—  dfSs  >  0,  i  =  l,...,n,  5  e  S1, 

(9) 

n 

X  —  Xj  >  0, 

i= 1 

(10) 

^Ps&s  >  Qu 

(11) 

seS 1 

-  Xi)  +  M( 1  -  Sst)  >  0, 

i  =  1, . . .  ,n,  s  e  S1,  t  e  S2, 

(12) 

n 

r-£y?>  o,  s  e  sl, 

i=l 

(13) 

^2  PstSst  >  Qi,  seS\ 

(14) 

teS2 


Xi,  Y  >  0,  x, ,  y*  €  Z,  Ss,  8st  €  {0, 1}, 

i  =  l,...,n,  s  e  Sl,  t  €  S2.  (15) 

The  binary  variables  Ss,  s  e  Sl,  and  Sst,  s  e  S\t  e  S2 
indicate  whether  a  certain  scenario  has  been  selected  to  be 
satisfied  in  day  1,  and  in  day  2  given  a  scenario  in  day  1, 
respectively.  Recall  that  the  requirements  in  day  2  must  be  sat¬ 
isfied  with  respect  to  each  possible  realization  of  the  demand 
scenario  in  day  1.  The  variables  yf,  i  =  l,...,n,  are  the 
recourse  variables:  the  flow  of  supply  to  (or  from)  D  Pi ,  given 
that  scenario  s  has  been  realized  in  day  1.  Recall  that  the 
value  of  some  y?  may  be  negative,  in  which  case  supply  is 
actually  taken  away  from  DP, .  The  constraints  in  (12)  deter¬ 
mine  the  demand  scenarios  in  day  2  to  be  satisfied,  where 
M  is  a  large  constant  compared  to  the  demand  data.  That 
is,  if  Sst  —  1  then  the  demand  of  scenario  t  is  to  be  sat¬ 
isfied  for  each  DPi,  given  scenario  .v  occurred  on  the  first 
day.  Since  inventories  cannot  be  negative,  we  require  that 
Xi  >  0,  i  =  1 Y  >  0.  The  following  example 
demonstrates  this  model. 
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Table  5.  Scenarios  in  day  1.  Table  7.  Conditional  transition  probabilities. 


Scenario 

DPi 

dp2 

DPi 

Pr. 

1 

10 

20 

30 

0.3 

2 

30 

60 

20 

0.4 

3 

20 

50 

10 

0.3 

EXAMPLE:  There  are  n  =  3  DPs,  and  on  each  day 
of  the  two-days  operation  there  are  three  possible  demand 
scenarios,  as  shown  in  Tables  5-7. 

The  probabilities  of  the  scenarios  in  day  1  are  given  in 
Table  5,  and  the  conditional  transition  probabilities  from 
scenarios  in  day  1  to  scenarios  in  day  2  are  presented  in 
Table  7.  The  demand  scenarios  of  day  1  and  day  2  are  given 
in  Tables  5  and  6,  respectively.  The  costs  of  a  unit  supply  are 
Cx  —  5  and  Cy  =  6  for  a  DP  and  the  depot,  respectively. 
The  required  minimum  probability  thresholds  are  Q\  =  1 
and  Q2  =  0.7,  for  day  1  and  day  2,  respectively.  An  opti¬ 
mal  deployment  (there  may  be  multiple  optima)  is  x  1  =  50, 
x2  =  80,  X3  =  60  with  X  —  xi  —  190,  and  Y  —  0.  If 
scenario  1  is  realized  in  day  1 ,  then  all  three  scenarios  can  be 
satisfied  in  day  2  (with  probability  1);  if  scenario  2  is  real¬ 
ized  in  day  1,  then  scenarios  2  and  3  are  satisfied  in  day  2 
(with  probability  0.4  +  0.3  =  0.7  as  required);  if  scenario  3 
is  realized  in  day  1,  then  once  again  all  three  scenarios  can 
be  satisfied  in  day  2  (with  probability  1).  The  values  of  the 
recourse  variables  if  scenario  1  is  realized  are:  y\  —  —10, 
y\  —  —50  and  V3  =  20,  which  means  that  some  of  the 
excess  supplies  in  DPj ,  i  —  1,2  is  distributed  to  D  P2  •  Sim¬ 
ilarly,  y\  —  0,  y\  =  —10,  y\  —  10,  and  yj5  =  0,  y\  =  —20, 
y\  —  0.  Notice  that  the  depot  is  initially  empty  (Y  =  0). 
This  comes  with  no  surprise  because  it  is  cheaper  to  store 
supplies  in  the  DPs  than  in  the  depot  (Cx  <  Cy )•  Absent 
transportation  costs  in  the  second  period,  it  would  always 
be  better  to  allocate  all  the  needed  supplies  to  the  DPs.  If 
the  costs  are  reversed,  that  is,  Cx  =  6  and  Cy  =  5,  then 
an  optimal  solution  is  x\  —  30,  x2  =  60,  x2  =  30  with 
X  —  ^?=1  xt  =  120,  and  F  =  70.  It  is  shown  below  that  the 
problem  can  be  decomposed  into  a  series  of  MKPs  and  that 
the  optimal  deployment  is  independent  of  the  actual  values 
of  the  costs  and  depends  only  on  the  ordinal  relation  between 
Cx  and  Cy. 

Next  we  show  that  the  two-period  problems  (8)— (15)  can 
be  decomposed  into  two  separate  problems,  each  is  an  MKP. 


Table  6.  Scenarios  in  day  2. 


Scenario 

DPi 

dp2 

DPi 

1 

30 

10 

10 

2 

20 

10 

50 

3 

10 

10 

20 

Scenarios  in  day  1 
Scenarios  in  day  2 

1 

2 

3 

1 

0.7 

0.3 

0.2 

2 

0.2 

0.4 

0.7 

3 

0.1 

0.3 

0.1 

The  period  1  problem,  which  is  labeled  <t>  [  and  is  equivalent 
to  (4)— (7),  is 

min  X 
s.t. 

Xj  —  d°8s  >0,  i  —  1, . . .  ,n,  s  e  Sl, 

n 

x-y>>o. 

1=1 

£pA  >  Qu 

seS 1 

Xj  >  0,  Xi  eZ,Ss  e  {0, 1),  i  =  1, . . . ,  n,  s  e  S1. 

Note  that  since  d-  >  0  for  all  i  and  s  e  Sl,  the  non¬ 
negativity  requirements  on  the  xt  variables  are  redundant. 
For  any  given  feasible  solution  x  —  (x\, . . .  ,xn)  of  <f>i,  the 
period  2  problem,  <J>2(x),  is, 

min  y(x) 
s.t. 

ysiix)  +  M(  1  -  8„)  >  «  +  d\  -  Xi), 

i  =  1, . . .  ,n,  s  e  51,  t  e  S 2, 

n 

y(x)  -  ^2  y'iW  >0,  jeS1, 

i=  1 

>  Qi,  seS1, 

reS 2 

y,s(x)  eZ,  8st  e  {0,1},  i  =  1, . . .  ,n,  s  e  S1,  t  e  S2. 

Note  that  y(x)  represents  the  total  sum  of  shipments  to  and 
from  the  depot  and  can  be  negative. 

Theorem  5.1  below  shows  that  an  optimal  solution  of  the 
two-period  problems  (8)— (15)  is  obtained  by  solving  0 1  and 
<l>2(x)  sequentially.  Recall  that  the  optimal  deployment  is 
independent  of  the  actual  values  of  Cx  and  Cy.  Moreover,  the 
optimal  total  amount  of  supply  (X+Y)  is  constant  for  a  partic¬ 
ular  problem  instance  and  is  independent  of  the  costs  Cx  and 
Cy.  This  property  suggests  the  model  is  robust  in  the  sense 
that  it  exempts  military  planners  from  specifying  exact  cost 
values  which  are  difficult  to  estimate.  They  need  to  express 
just  ordinal  preferences  between  two  logistic  options.  Also, 
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Proposition  5.2  below  indicates  that  the  integrality  constraints 
on  the  x  and  the  y  variables  in  these  problems  are  redundant 
if  the  demand  values  d  are  integer.  We  leave  the  proofs  of 
these  results  to  the  Appendix. 

THEOREM  5.1:  Let  x  be  an  optimal  solution  for  <t>i  and 
y(x)  be  the  corresponding  optimal  objective  value  of  O 2 ( -V ) , 
and  let  X  =  Y2'i=i  *i  anc*  f (*)  —  max{0,  v(a)}.  Then, 

(i)  if  Cx  >  Cy  then  (. X ,  Y (Jc ))  is  optimal  for  (8)— (15); 

(ii)  if  Cx  <  Cy  then  (X  +  Y(x),0)  is  optimal  for 
(8)— (15). 

For  each  of  the  problems  <t>i  and  02(a)  denote  by  MlP-<t>i 
and  MIP-TlIx),  respectively,  the  former  problems  where  the 
integrality  requirements  are  not  imposed  on  the  x  and  the  y 
variables. 

PROPOSITION  5.2:  Let  ds  and  d'  be  integral  demand  vec¬ 
tors  for  each  s  e  Sl  and  t  e  S2,  respectively.  Then,  there 
always  exist  optimal  solutions  of  MIP-<J>i  and  of  MIP-TMa) 
which  are  integral. 

We  have  already  argued  in  Section  3  that  <t>i  is  essentially 
an  MKP.  We  wish  to  show  now  that  a  solution  for  <J>2(x)  is 
obtained  by  solving  |  S'1 1  MKPs.  For  a  given  period  1  solution 
x  and  a  given  period  1  scenario  s  e  Sl,  consider  problem 
d>*(x)  below: 

n 

min  ^y-  (x) 

1  =  1 
s.t. 

y'l  (x )  +  M(  1  -  Ssr)  >  (df  +d‘-xi),  i  =  l,...,n,te  S2, 
Pstsst  >  Qi, 

teS2 

y-(x)  e  Z,  8st  e  {0, 1},  i  =  1, . . .  ,n,  t  e  S2. 

Clearly,  each  one  of  the  .S’1 1  problems  <t>j(x)  is  an  MKP. 

PROPOSITION  5.3:  For  a  given  period  1  solution  x,  let 
ys(x )  =  (yj(x), . . . ,  y'n(x))  be  an  optimal  solution  of  <J>,(x), 
5  e  S1,  and  let  s  e  argmax{^"=1  y-(x)  :  s  e  51}.  Then, 
y(x),  an  optimal  solution  of  <62 (a)  is  obtained  by, 

n 

y(x)  =  (a). 

i=i 

PROOF:  The  objective  in  problem  is  to  select  the 

minimum  y(x)  value  that  satisfies  the  probability  threshold 
requirement  (14)  for  all  period  1  scenarios  s  e  Sl.  This 
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minimum  is  attained  if  for  each  problem  OjCx),  s  e  Sl, 
1  (a)  is  minimized.  y(x)  is  the  maximum  among  these 

values.  □ 


6.  SUMMARY  AND  CONCLUSIONS 

In  this  paper  we  present  a  new  knapsack-related  combina¬ 
torial  problem  termed  minmax  multidimensional  knapsack 
problem  (MKP)  that  is  motivated  by  a  military  logistics 
problem.  The  logistics  problem  is  to  determine  an  optimal 
deployment  of  inventories  that  satisfies  certain  operational 
requirements  in  a  two-days  scenario.  Recourse  opportunities, 
which  are  crucial  in  combat  related  military  environments, 
are  explicitly  incorporated  in  our  model.  We  show  that  the 
resulting  two-period  stochastic -programming  problem  can  be 
solved  by  solving  a  series  of  MKPs.  A  practically  efficient 
algorithm  is  developed  for  solving  the  MKP. 

The  above  result  can  be  generalized  to  the  multi-period 
problem  with  recourse.  The  multi-period  problem  can  be 
decomposed  into  T  separate  problems,  <t>i,  <t>2, . . . ,  ‘Fr,  that 
correspond  to  periods  1,2 , ,T ,  respectively.  We  observe 
that  to  solve  the  multi-period  problem  we  need  to  opti¬ 
mize  only  the  first  and  the  last  periods’  allocations.  For  all 
intermediate  periods  we  need  to  find  just  feasible  solutions 
that  satisfy  the  corresponding  responsiveness  chance  con¬ 
straints.  This  observation  follows  from  the  replenishment 
policy,  where  all  tapped  safety  stocks  are  replenished  in 
the  subsequent  periods,  and  the  assumption  regarding  free 
transshipment. 

Possible  extensions  of  the  logistics  model  involve  addi¬ 
tional  constraints  such  as  zero  safety  stock,  which  excludes 
back-orders,  and  no  transshipment,  which  implies  that  the 
recourse  variables  must  be  non-negative.  Another  possible 
extension  may  be  to  incorporate  possible  (stochastic)  inter¬ 
diction  on  the  transportation  of  supplies  from  depot  and 
among  the  DPs.  These  extensions  may  generate  extensions 
to  the  MKP  where  side  constraints  are  incorporated  in  the 
knapsack-like  setting. 


APPENDIX 

We  prove  now  Theorem  5.1  and  Proposition  5.2.  First  we  show  that  the 
resulting  sequential  optimal  solutions  for  O  i  and  $2  provide  an  optimal  solu¬ 
tion  for  the  combined  problem.  To  prove  this  result  we  need  the  following 
lemma. 

LEMMA  Al:  Let  x'  and  x"  be  two  feasible  solutions  for  Oi  such  that 
x"  =  x'i  +  of,- ,  a,-  e  R,  i  =  1 Let  y(V)  and  y(x")  be  the  cor¬ 
responding  optimal  objective  values  of  an^  respectively. 

Then,  y(x ")  =  y(x ')  -  £"=1  a,-. 

PROOF:  We  start  by  showing  that  y  {x")  >  y  (V)  —  Y^i=i  ai  •  Let  x  be  any 
feasible  solution  of  <l>i  and  let  a*’1  (x)  =  d ?  +  d*  —  x;  for  any  i  =  ,n, 
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s  6  S1,  t  e  S2.  Then,  for  each  i  €  S1  we  define 

As{x)  =  {as-'(x)  =  (af-' :/eS!|. 

Let  as(x )  be  an  optimal  solution  to  the  MKP  on  As(x)  and  observe  that 
y(x)  =  max{fls(x)  :  s  e  S'1}-  Then  it  suffices  to  show  that  for  each 
s  e  S1,  the  inequality  as(x")  >  as{x ')  —  E?=i  ai  holds.  This  is  because 
y(x")  >  as(x ")  >  as(x ')  -  E?=i  «i  for  each  s  e  s'  and  = 
max{fl,(jc/)  :  s  e  S1}.  Let  s  be  any  scenario  in  S’1.  Assume  on  the  contrary 

that  as(x")  <  as(x ')  -  E?=i  then  E7=i («?(*")  +  «i)  <  £?=i  «?(*')• 

We  show  below  that  as(x")  +  a,  where  a  =  (ar,*)"=1 ,  is  a  feasible  solution 
for  the  MKP  on  As(x'),  contradicting  the  optimality  of  as{xr). 

Let  as( x")  be  an  optimal  solution,  and  hence  a  feasible  solution,  for  the 
MKP  on  As  (x").  Therefore,  there  exists  a  feasible  subset  W  c  S2  such  that 
for  each  t  e  W  and  for  each  i,  the  inequality  a f  (x")  >  a^(x")  holds.  Also, 
for  each  t  e  W  and  for  each  i,  by  the  definition  of  (x)  and  of  x"  the 
equality  <2{-  (x")  =  (xr)  —  at  holds.  Hence,  for  each  t  e  W  and  each  i, 
af(x")  >  a5/ ( x")  =  a/ix')  -  implying  a?  (x")  +  >  a5/ (xf).  Thus, 

as  (x,;)  +  a  is  a  feasible  solution  for  the  MKP  on  the  set  As  (V),  of  value  less 
than  as(x'),  contradicting  the  optimality  of  as(x'). 

Showing  the  other  direction  that  y(x")  <  y(x')  —  E"=i  ai  can  done 
in  the  same  way  by  showing  that  y(x')  >  y(^/r) —  E?=  1  (— ai)-  Thus,  by  the 
two  directions  we  get  that  y(V0  =  y(V)  —  E?=i  ai  as  required.  □ 

The  next  theorem  states  that  the  sequential  procedure  for  obtaining  the 
logistics  deployments  in  periods  1  and  2  is  globally  optimal. 

THEOREM  5.1:  Let  x  be  an  optimal  solution  for  <J>i  and  y(x)  be  the 
corresponding  optimal  objective  value  of  <I>2(a),  and  let  X  =  E?=i  and 
Y (x)  =  max{0,  y(^)}-  Then, 

(i)  if  Cx  >  Cy  then  (X,  F(^))  is  optimal  for  (8)— (15); 

(ii)  if  Cx  <  Cy  then  ( X  +  T(x),0)  is  optimal  for  (8)— (15). 

PROOF:  (i)  To  show  the  correctness  of  the  theorem  in  the  case  of 
Cx  >  Cy  it  suffices  to  show  that  for  any  feasible  solution  x  of  such  that 
X  =  X  +  E"=i  ai  with  E/Li  ai  —  where  X  =  E”=i  xi->  the  inequality 
CxX  +  Cy Y (x)  >  CxX  +  CyY(x)  holds. 

Suppose  T(x)  =  0,  which  yields  y(x)  <  0.  If  F(Jt)  =  0  then  the  above 
inequality  follows  from  the  fact  that  X  >  X.  Otherwise,  Y (x)  >  0,  implies 
y(L)  >  0.  By  Lemma  AI,  y(L)  —  E?=i  ai  ~  y(x)  —  which  implies 
Ya=\  “i  i  Kx)  =  Y(x).  Then, 


CxX  +  CyY(x )  =  CxX  +  Cx  X  a‘  +  Cy  ■  0 
i= 1 

>  CxX  +  CXY(X)  >  CxX  +  CyY (x). 

Suppose  now  Y  (x)  >  0,  which  implies  that  Y  ( x )  =  y(x)  >  0.  By  Lemma 
AI,  since  X  >  X,  we  get  Y (x)  =  Y (x)  —  ^2"=1  o', .  Then, 


CxX  +  CyY(x)  =  CxX  +  Cx  J2  “i  +  CyY(x)  -  CyJ^  «< 

i=l  i=l 

n 

=  CxX  +  CyY(x)  +  (CX  -  CY)J2a-  >  CxX  +  CYf(x). 

1  =  1 


(ii)  In  the  case  of  Cx  <  Cy  we  only  need  to  show  that  CxX  +  Cy  F(x)  > 
Cx  (X  +  F (x))  for  any  feasible  solution  x  of  <3>i .  The  above  follows  directly 
from  Cx  <  Cy .  □ 

We  next  turn  to  show  that  the  integrality  requirements  on  the  x  and  the 
ys  variables  in  our  IP  formulations  are  redundant.  For  each  of  the  problems 
C>1  and  <E>2(x)  denote  by  MIP-<Li  and  MIP-<I>2(x),  respectively,  the  former 
problems  where  the  integrality  requirements  are  not  imposed  on  the  x  and 
the  ys  variables. 

PROPOSITION  5.2:  Let  ds  and  d 1  be  integral  demand  vectors  for  each 
s  e  S’1  and  t  e  S2.  Then,  there  always  exist  optimal  solutions  of  MIP-<3>i 
and  of  MIP-<3>2Cr)  which  are  integral. 

PROOF:  The  proposition  follows  immediately  from  the  fact  that  the 
8t  variables  are  {0,  l}-variables  and  the  integrality  of  the  demand  vectors. 
Observe  that  in  problem  MIP-<t>i  we  seek  for  a  subset  of  scenarios  (this  is 
since  8S  is  a  binary  variable)  for  which  the  required  integral  demands  will  be 
satisfied.  The  above  coupled  with  the  minimality  of  X  imply  the  integrality 
of  the  x’s.  Similarly,  it  can  be  shown  that  MIP-d>2(x)  is  integral  for  any 
integral  vector  x  and  integral  demand  vector.  □ 
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